第 1 卷 第 1 期 原子 核 物 理 评 论 Vol. 1, No. 1 
2024 年 2 月 Nuclear Physics Review Feb. , 2024 
文章 编 1007-4627(2024) 01-0001-07 


CNPC2023+ 多 重 物 理 约束 下 Duflo-Zuker 核 质量 模型 的 研究 


叶 伟 虎 !， 钱 以 斌 ! 
(1. 南京 理工 大 学 物理 学 院 ， 南 京 210094) 

摘要 : 原子 核 质量 在 核 物 理 和 核 天 体 研究 中 扮演 着 十 分 重要 的 角色 。 尽 管 当前 的 理论 模型 可 较 好 的 描述 已 知 
核 素 的 质量 ， 但 不 同 模型 在 外 推 上 仍然 存在 许多 和 争议， 尤其 是 丰 中 子 区 的 预言 。 本 文 回顾 了 我 们 近年 来 在 核 
质量 模型 外 推 研究 的 主要 成 果 ， 即 使 用 严格 的 多 目标 优化 算法 将 核 素 质量 差 a 衰变 能 和 局 域 质 量 关 系 作为 限 
制 条 件 ， 缓 和 质量 模型 中 的 过 拟 合 问题 ， 在 一 定 程度 上 提升 模型 的 外 推 能 力 。 
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1 引言 邻 核 素 间 的 质量 关系 间接 获得 未 知 核 素质 量 ， 其 理论 靖 
Lk bb 度 约 在 150keV 以 内 ， 包 括 有 著名 的 Garkev-Kelson(GK) 
eae 
VH. NE. BRONTE. Gewuy guo SEARO 等 等 。 国 内 课题 组 在 局 域 质量 关系 中 
REO. AODH SAME am 。 修 出 了 许多 重要 的 发 现 和 拓展 ， 包 括 对 GK 关系 做 出 进 
格 计算 的 重要 输入 ,这 对 理解 21 世纪 尚未 解决 的 11 一 步 推广 并 确定 最 优化 外 推 路 径 的 研究 ARA 
个 重大 物理 难题 之 三 有 着 重要 积极 意义 ， 即 “ 字 宙 中 从 ”外 推 公式 的 研究 。 “等 等 。 关 于 局 域 硕 全 关系 的 更 多 
铁 到 钠 的 元 素 是 如 何 产生 的 ” 因此, 质量 测量 一 直 是 核 细节 可 查阅 最 新 综述 “。 虽 然 宏 观 -微观 模型 的 物理 思 
科学 研究 的 热点 。 想 广 度 和 深度 要 lin 但 在 一 定 质 量 数 区 域 
实验 上 可 通过 直接 或 间接 测量 的 方法 获得 约 2500 ”站 局 域 型 的 结合 能 精度 远 远 高 全 局 型 。 近 年 来 ， 随 着 
个 核 来 的 抽 量 ， 市 不 同 理论 预计 给 出 在 折 然 界 由 可 在 在。 核 物理 和 机 器 学 习 交叉 学 科 的 兴起 ， 核 质量 模型 的 精度 
的 核 素数 目 约 是 7000-10000Kt。 目 前 由 于 这 些 未 知 核 。 获得 巨大 改善 S, IERI 100keV 以 内 ， 即 满足 核 合 
素 集中 在 质子 滴 线 区 、 丰 中 子 区 和 超重 核 区， 目标 核 寿 。“ 半 站 所 二 的 精度 需求 。 作 是 ， 近 几 年 间 AANE 
命 短 、 产 生 截面 小 ， 精 确 测量 其 质量 存在 着 许多 技术 上 出 了 一 个 尖锐 的 问题 模型 在 忆 知 区 域 核 过 的 精度 和 未 
的 困难 。 在 可 预见 的 将 来 ， 不 稳定 核 素质 量 仍 需 依赖 理 VORNE ENERE? 
iex S. 不 幸 的 是 ，Sobiczewsk 等 人 基于 当前 流行 的 数 十 个 
原子 核 质量 的 理论 模型 通常 可 分 为 两 类 。 一 类 是 全 。 质量 模型 的 系统 性 研究 时 发 现 模型 在 已 知 区 域 核 素 的 精 
局 型， 另 一 类 为 局 域 型 。 全 局 型 中 精度 较 好 的 是 宏观 -。 度 和 未 知 区 域 的 预言 能 力 的 相关 性 较 小 BT。 这 在 一 定 
微观 和 纯 微观 模型 ， 宏 观 -微观 模型 的 宏观 项 指 的 是 描 “程度 上 表明 单一 追求 结合 能 的 精度 可 能 并 不 会 提升 模型 
述 结合 能 的 平滑 项 ， 微 观 一 般 指向 壳 层 和 对 能 修正 。 纯 。 计算 未 知 核 的 预言 能 力 .另外 ,质量 模型 还 存在 着 参数 化 
微观 模型 的 经 典 思 想 为 给 定 核 子 相互 作用 势 场 ， 在 非 。 模型 的 内 部 矛盾 , 即 过 拟 合 现象 。 比 如 , 即使 DZ-type Jf 
相对 论 下 为 求解 苹 定 询 方 程 获得 结果 ， 在 相对 论 下 则 。 量 模 型 计算 结合 能 的 精度 较 高 ， 但 质量 相关 量 a 衰变 能 
自然 为 求解 狄 拉克 方程 。 较 为 常用 的 全 局 型 包括 Duflo- 的 精度 相对 较 低 B9。 不 仅 如 此 ， 不 同 质量 模型 的 外 推 关 
es 1 模型 、 有 限 力 程 小 液 滴 模 型 中 和 平均 场 理 。” 异 显著 ， 在 中 子 滴 线 附近 差异 能 达到 几 十 个 MeVB439 。 
习 等， 其 标准 差 约 在 300-600keV 人 I。 局 域 型 则 利用 相 。 ”这 些 现象 在 一 定 程度 上 表明 了 质量 模型 存在 不 可 忽视 的 
过 拟 合 问题 。 虽 然 过 拟 合 的 模型 能 够 完美 重 现 拟 合 样本 
omn. BRERA OME UNUD, PARREREN 集 的 数据 ， 但 是 在 预测 样本 集 外 的 数据 时 ， 误 差 将 随 着 
mre mum (30922010312) _ 外 推 距离 的 增 大 而 增 大 。 
ene n ts a 通常 ， 缓 和 过 拟 合 的 主要 方案 是 提供 更 多 的 拟 合 样 
UL ETHER 本 ,如 此 可 避免 算法 在 某 一 区 间 内 过 度 “ 学 习 ” 进而 达 
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到 自由 参数 全 局 最 优 而 非 局 部 最 优 。 因 此 ， 我 们 在 核 质 
量 领域 中 引入 多 目标 优化 算法 ， 将 a 衰变 能 和 GK 关系 
作为 两 个 重要 条 件 去 约束 质量 模型 。 本 文 肯 在 通过 引入 
多 个 物理 约束 条 件 ， 提 升 质量 模型 在 远 距 离 外 推 时 的 预 
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2 模型 介绍 和 算法 


Duflo-Zuker (DZ) 系列 是 目前 研究 中 使 用 率 较 高 的 
质量 模型 , 原因 在 于 其 具备 较 高 的 理论 精度 和 外 推 性 能 ， 
其 计算 成 本 也 要 远 远 的 小 于 其 他 同等 精度 的 质量 模型 。 
DZ 系列 模型 包括 DZ10、DZ28 和 DZ33 等 B738， 本 文 
采用 DZ10 模型 为 例 ， 给 出 其 多 目标 优化 下 的 结果 。 


24 DzZ10 模型 


DZ 模型 由 宏观 部 分 的 液 滴 项 和 微观 壳 模 型 
monopole 哈密 顿 项 组 成 。DZ10 模型 对 一 个 原子 核 的 基 
态 结合 能 可 以 如 下 表示 : 


M 
B Epzio =a; Vç + a (M + S) 一 043 — 十 a4 Vp Ea as Vy 
p 


S 
+ agVprg + 4584 一 ag > + d954 + d49d4 (1) 
p 


其 中 : a, 是 由 拟 合 确定 的 自由 参数 ，aj、a4、a5 和 ag 分 
别 代表 液 滴 模 型 里 的 库仑 能 、 对 能 、 对 称 能 和 对 称 能 修 
正 能 ， 其 余 项 是 微观 价 核子 修正 。 


2.2 多 目标 优化 算法 


人 们 在 职场 和 生活 中 往往 会 面临 许 许多 多 的 优化 问 
题 。 比 如 ， 公 司 收入 利润 最 大 化 、 生 活 成 本 最 小 化 等 问 
题 。 然 而 ， 正 如 事物 之 间 具 有 相关 性 ， 对 某 一 目标 要 求 
收益 最 大 化 则 必然 会 导致 其 他 相关 目标 收益 降低 ， 因 此 
需要 对 不 同 目标 间 的 冲突 关系 做 出 阐述 ， 这 类 同时 考虑 
两 个 目标 以 上 的 优化 问题 称 为 多 目标 优化 。 多 目标 优化 
算法 被 广泛 应 用 在 经 济 学 、 工 程 学 等 学 科 B?41。 多 目标 
优化 算法 定义 了 帕 累 托 前沿 解 (Pareto ffont)， 又 被 称 为 
非 劣 解 解 集 。 非 劣 解 解 集中 包含 了 一 系列 基于 不 同 目标 
补偿 关系 所 确定 的 解 ， 非 劣 解 的 任何 一 个 解 儿 乎 不 可 能 
使 所 有 目标 同时 达到 最 优 的 效果 ， 其 核心 思想 是 从 各 个 
冲突 的 目标 函数 条 件 下 寻找 尽 可 能 多 的 补偿 解 ， 最 终 提 
供用 户 尽 可 能 多 的 选择 。 不 失 一 般 性 ， 一 个 多 目标 优化 
问题 可 以 如 下 表示 : 


min F(x) = (fi), .... f Q0) 
s.t. xcaQ 


lb € x; < ub. (2) 


其 中 : f(x) 是 目标 函数 ; x = (x1,xX2,.… Xn) 为 决策 变 
量 ;每 个 决策 变量 被 限制 在 上 界 和 下 界 的 解 空间 中 。 由 
于 多 目标 优化 中 涉及 多 个 目标 函数 的 比较 ， 因 此 无 法 像 
单 目标 优化 中 单独 比较 一 个 函数 值 就 可 确定 个 体 的 优 
劣 。 帕 累 托 支配 关系 是 多 目标 优化 算法 的 核心 要 义 ， 殿 
定义 解 空间 中 个 体 的 优 劣 关系 。 假 设 解 空间 中 有 两 个 个 
fk s Flt, 个 体 支配 ( 优 于 ) Ma 1, FAM 4X Vi = 
1.2, ...,m,f;(s) € OE 3je[1,2, ... m], HE fs) « 
f; Me 


3 结果 及 分 析 


3.1 a 衰变 能 的 约束 
前 ，non-dominated sorting genetic algorithm-II 
(NSGA-II) 仍然 是 第 二 代 进 化 算法 中 使 用 率 最 大 、 引 用 
率 最 高 的 算法 此 人。 因此 ， 本 文选 择 此 算法 对 结合 能 和 
a 衰变 能 进行 多 目标 优化 ,a 衰变 能 可 由 母 核 和 子 核 的 质 
BEAM, BO, = BE(Z-2, N C2) BE, - BE(Z, N). 
目标 函数 的 构造 由 理论 值 和 实验 值 的 标准 差 给 出 : 


Zi (ken = Ktheo.)2 
Wi =EN mnm (3) 
N 


其 中 : o, 为 计算 理论 值 和 实验 值 的 标准 差 公 式 ， 六 为 
实验 样本 数 ， 结 合 能 数量 为 2340 个 ， 偶 偶 核 a 衰变 能 
数量 为 275 个 ; k 是 对 结合 能 或 其 他 物理 量 的 代称 。 结 
合 能 和 a 衰变 能 的 实验 数据 取 自 国内 的 原子 核 质 量 评估 
AME2020[4。 图 1 为 通过 NSGA-II 算法 得 出 的 DZ10 模 
型 非 劣 解 结果 图 。 横 、 纵 轴 分 别 代 表 结 合 能 和 a 衰变 能 
的 标准 差 。 首 先 ， 图 1 中 的 左上 第 一 个 解 可 理解 成 既是 
结合 能 精度 最 优 解 又 是 a 衰变 能 精度 最 劣 解 ， 而 右 下 的 
最 后 一 个 解 则 相反 。 我 们 可 知 质量 模型 中 结合 能 和 a E 
变 能 的 关系 是 相互 冲突 的 ， 即 结合 能 精度 的 提升 是 以 劣 
化 a 衰变 能 精度 为 代价 。 这 告诉 我 们 以 往 的 单 目 标 算 法 
在 拟 合 结合 能 确定 参数 时 ， 是 以 质量 为 最 高 优先 级 的 局 
部 最 优化 ， 而 不 是 考虑 结合 能 和 o 衰变 能 多 个 物理 量 的 
全 局 最 优 ， 这 可 能 会 使 得 过 拟 合 问题 加 重 ， 劣 化 模型 在 
拟 合 训 练 集 之 外 的 预言 。 其 次 ， 在 考虑 了 a 衰变 能 的 数 
据 集 后 ， 可 获得 一 系列 计算 a 衰变 能 更 占 优 的 解 。 实 际 
E, a 衰变 能 是 a 衰变 寿命 的 一 个 重要 决定 因素 , 以 往 的 
计算 中 发 现 a 衰变 能 0.01MeV 的 不 确定 度 最 终 可 影响 a 
半衰期 一 个 数量 级 ， 这 些 依赖 模型 本 身 获得 的 a 衰变 能 
占 优 解 可 以 对 预测 重 核 和 超重 核 的 a 衰变 能 提供 一 定 支 
持 。 近 年 来 ， 由 于 不 同 理论 模型 在 预言 实验 值 之 外 的 区 
域 时 往往 差异 过 大 ， 因 此 理论 学 家 则 需要 给 出 模型 的 系 
统 和 统计 不 确定 度 来 对 模型 自身 的 不 确定 度 给 出 定量 描 
述 。 基 于 这 几 十 个 非 劣 解 ， 我 们 可 从 给 出 模型 的 统计 不 
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表 1 DZ10 模型 的 统计 不 确定 度 和 质量 最 优 解 参 数 ， 单 位 


0.44 © DZ10 MeV. 
0.43 
S 0.42 
[^] 
= on 参数 I II 
3 
a 
AN a, 0.705+0.001 — 0.705 
0.39 a,  17735x008 17747 
0.38 a, 16.207 土 0.093 16251 
0.0 0.5 1.0 1.5 
ope(MeV) ay 5.906+0.58 6.102 
a, — 37.303+0.097 37.356 
图 1 (在 线 彩 图 ) 使 用 NSGA-II 算法 优化 DZ10 质量 模型 的 非 ag 52.757 + 0.568 — 52.661 
劣 解 结果 图 a, 0.438 土 0.009 0.463 
ay — 2.006-0.041 — 2.104 
MER. R 1 为 DZ10 模型 参数 的 不 确定 度 。 第 二 列 为 参 a,  0.022+0.0005 0.021 
数 的 统计 不 确定 度 ; 第 三 列 为 质量 最 优 解 的 参数 旨 在 给 aio 39.602 £ 0.591 41.48 
出 对 比 oo, 0.397 + 0.011 0.439 
- "m 2 、 " o 0.66 + 0.083 0.562 
基于 非 劣 解 解 集 的 统计 性 ， 图 2 系统 性 的 给 出 了 新 = 


Ww 


致 的 讨论 。 


的 核 素 边界 、 双 中 子 分 离 能 和 非 劣 解 的 外 推 差 异 等 物理 。 条 件 ， 对 质量 模型 做 出 更 
信息 。 第 一 , 液 滴 线 可 利用 核子 分 离 能 小 于 零 确定 , 由 于 
偶偶 核 往往 比 奇 A 核 稳 定 ， 所 以 通过 双核 子 分 离 能 确定 

的 液 滴 线 包含 了 通过 单 核 子 分 离 能 确定 的 液 滴 线 。 图 中 ， AMex = M(N +2,Z-2)- M(N,Z)+ M(N.Z - 


P(S,, < 0) 为 通过 双 中 子 分 离 能 确定 的 该 核 素 是 否 为 非 —M(N +1,Z —2)+ MIN -L,Z)- MIN -2,Z- 1) 
束缚 核 的 概率 ，P(S,p < 0) 为 通过 双 质 子 分 离 能 确定 的 (4) 


该 核 素 是 否 为 非 束 缚 核 的 概率 。 第 二 ， 我 们 系统 的 计算 
了 非 劣 解 的 双 中 子 分 离 能 ， 并 取 加 权 平 均 ， 大 小 通过 图 
左 侧 数值 颜色 棒 表 示 。 第 三 ， 我 们 采用 非 劣 解 解 集中 
的 质量 最 优 解 和 次 质量 最 优 解 对 应 参数 组 ， 即 图 1 中 最 
上 方 相 邻 的 两 个 解 ， 并 计算 其 之 间 的 的 结合 能 差异 ， 在 
右 下 子 图 表示 。 此 子 图 清晰 表明 在 有 实验 值 和 丰 质 子 区 
时 ， 外 推 的 差异 不 超过 0.5MeV， 但 是 在 超重 核 区 和 让 
子 区 时 ， 差 异 则 会 明显 增 大 。 我 们 考虑 到 这 两 个 解 具 
备 明显 偏向 质量 精度 的 特征 ， 因 此 可 能 会 导致 大 规模 外 
E 时 过 拟 合 现 象 加 重 ， 最终 劣 化 模型 外 推 的 差异 。 于 是 ， 
我 们 从 图 1 的 中 间 区 域 任 取 两 个 解 计算 外 推 的 差异 。 通 
过 比 对 图 2 中 子 图 和 图 3， 可 知 中 间 区 域 的 外 推 差异 要 


其 中 : AMgx 是 相 邻 6 个 核 素 加 加 减 减 抵消 后 的 不 确定 
常数 , 通常 在 +150keV 以 内 。 我们 基于 式 4 从 AME2020 
中 提取 了 几 百 组 AMEE 的 实验 值 ， 并 将 理论 DZ10 BE 
型 得 到 的 AM he 与 其 做 标准 差 。 实 际 中 ， 基 于 剩余 核 
子 相 互 作 用 相互 抵消 假设 ， 衍 生出 了 许多 版 本 的 GK 关 
系 ， 可 多 达 数 十 个 。 本 文 简单 采用 了 式 4 将 其 作为 可 实 
现 案例 ， 表 明了 利用 GK 关系 多 目标 限制 质量 模型 的 可 
行 性 ， 当 然 也 可 将 式 4 替换 成 其 他 版 本 的 GK 关系 ， 这 
还 有 待 未 来 做 出 进一步 的 检验 。 图 4 为 多 目标 优化 算法 
下 以 结合 能 、a 衰变 能 和 GK 关系 给 出 的 非 劣 解 解 集 图 。 
青色 方 框 包 庄 的 为 质量 最 优 解 。 上 文 已 经 谈论 过 多 个 物 


理 量 的 平衡 解 在 外 推 上 可 能 更 可 靠 ， 但 此 说 明 还 仍 待 进 
明显 的 小 于 质量 优 解 组 世 ， 表 明 利 用 平衡 结合 
din iHd aun NAR UR isla 一 步 确认 ， 下 文 将 对 此 关键 问题 做 出 细致 分 析 。 面 对 如 


We 比 罗 的 非 和 解 ， 如 何 从 中 寻找 出 更 有 物理 让 义 、 更 阴 给 


3.2 Garvey-Kelson 关系 的 约束 出 可 靠 外 推 值 的 解 是 十 分 关键 的 问题 。 由 于 丰 中 子 核 素 
在 前 文 讨论 的 基础 上 , 可 以 对 质量 模型 的 约束 条 件 。 其 备 较 大 的 同位 旋 不 对 称 度 ， 对 理论 和 实验 而 言 预测 和 
做 进一步 的 探索 , 即 通 过 引入 局 域 质量 公式 来 约束 模型 。 ”测量 这 些 核 素 质量 都 是 巨大 的 挑战 ， 且 不 同 理论 模型 的 
本 文 将 Garvey-Kelson 关系 作为 限制 核 质量 模型 的 另外 AREER PFEP TKR. 

一 个 条 件 。 过去, J. Barea 等 人 就 提出 可 对 全 局 型 模型 拟 首先 ， 我 们 从 近 几 年 的 实验 中 提取 出 尚未 包括 在 
n rem e Gan AME2020 "HAUSE FAY iE us OOS), PR ER DZ10 
邻 核 的 连贯 性 和 外 推 时 的 预测 性 四 1, 但 文中 讨论 的 实践 。 ”的 质量 最 优 解 计算 这 些 核 的 标准 差 并 将 其 作为 一 个 理论 
方法 较 粗 糙 且 元 余 。 之 后 ， 全 局 型 模型 往往 将 GK 关系 ” ”基准 值 ， 若 其 余 的 非 劣 解 给 出 比 此 基准 值 精度 更 好 的 标 
作为 后 验 样本 ， 而 非 限制 条 件 。 借 助 多 目标 优化 算法 的 佳 差 时 ， 将 保存 这 些 非 劣 解 并 归 入 到 潜力 解 解 集 。 事 实 
优势 ， 我 们 将 Garvey-Kelson 关系 作为 一 个 严格 的 限制 Eo 非 劣 解 的 数目 高 达 几 百 个 , 但 筛选 出 来 的 只 有 几 十 个 
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图 2 (在 线 彩 图 ) DZ10 模型 非 劣 解 解 集 的 统计 结果 图 。 上 、 下 两 侧 分 别 是 以 概率 形式 给 出 的 双 质 子 、 双 中 子 滴 线 ， 绿 色 方 块 代 


K AME2020 中 有 质量 实验 数据 的 核 素 ;中 间 的 折线 为 p 稳定 线 ， 颜 色 圆 圈 为 非 劣 解 下 双 中 子 分 离 能 的 平均 值 ; 右 下 角子 


ABE (MeV) 
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图 3 (在线 彩 图 ) 左 子 图 (a) 为 从 图 1 的 中 间 区 域 的 两 个 解 对 
应 参数 组 所 给 出 的 结合 能 差异 。 右 子 图 b) 为 从 图 4 的 
中 间 区 域 的 两 个 解 对 应 参数 组 所 给 出 的 结合 能 差异 。 
潜力 解 , 这 些 潜力 解 自 然 可 以 提供 一 些 有 用 信息 。 表 2 从 
局 部 和 整体 两 方 对 潜力 解 给 出 系统 性 分 析 结 果 ， 第 二 至 
第 四 行 针 对 局 部 丰 中 子 核 数 据 ， 最 后 三 行 针 对 整体 拟 合 
样本 数据 。 其 中 : 第 二 行为 质量 最 优 解 计算 得 出 的 标准 
差 ， 第 三 行为 潜力 解 解 集中 所 给 出 最 好 的 标准 差 ， 第 四 
行为 所 有 潜力 解 结合 能 标准 差 的 统计 分 析 ， 第 一 项 为 标 
准 差 的 平均 值 , 第 二 项 为 标准 差 秒 速 浮动 情况 ; 第 五 、 六 、 
和 七 行 分 别 为 潜力 解 解 集 计算 拟 合 样 本 中 结合 能 、a E 
变 能 和 GK 关系 的 统计 人 分析。 首先， 最 好 潜力 解 计算 的 
丰 中 子 核 质量 标准 差 远 远 高 于 质量 最 优 解 结果 ， 约 提升 
70% 的 精度 。 其 次 ， 观 察 表 中 第 四 行 可 知 ， 这 些 数 据点 
的 变化 浮动 较 小 ， 且 潜力 解 解 集 对 a 衰变 能 和 GK 关系 
目标 的 平均 值 集 中 于 各 自 的 中 间 区 域 ， 即 平衡 区 域 。 这 
在 一 定 程度 上 验证 了 上 文 的 论点 ， 平 衡 多 个 物理 量 下 的 
全 局 解 更 有 可 能 给 出 稳定 可 靠 的 外 推 值 。 


上 文 在 考虑 结合 能 和 a 衰变 能 的 关系 下 ， 采 用 了 
平衡 两 者 关系 的 两 个 全 局 最 优 解 计算 其 之 间 的 结合 能 差 


图 是 DZ10 基于 非 劣 解 解 集中 质量 最 优 解 和 次 质量 最 优 解 所 对 应 参数 组 的 结合 能 差异 图 。 


表 2 潜力 解 解 集 的 统计 分 析 结 果 。 二 至 四 行为 潜力 解 解 集 
对 拟 合 训练 集 之 外 的 丰 中 子 核 结 果 ; 后 三 行 分 别 代表 潜 
力 解 解 集 对 训练 集 内 的 结合 能 、a 衰变 能 和 GK 关系 下 
的 结果 ， 单 位 MeV。 


物理 量 数值 

oe 0.435 

one 0.091 

ofS 0.304 + 0.095 

Opp 0.823 + 0.244 

og, 0.362 + 0.016 

cju, 0-066 + 0.001 
异 ， 发 现 结合 能 在 外 推 时 的 差异 已 经 减 小 ， 在 某 种 程度 
上 表明 过 拟 合 得 到 一 定 改善 。 自 然 ， 我 们 需要 考虑 在 加 
入 了 结合 能 、a 衰变 能 和 GK 关系 下 ， 即 更 多 物理 数据 


的 约束 下 ， 提 取 的 非 劣 解 的 结合 能 差异 情况 。 图 3 的 子 
图 (b) 为 从 潜力 解 解 集 ( 即 中 间 区 域 ) 任 挑 两 组 解 计算 
的 结合 能 差异 图 。 相 比 于 子 图 a), FA Coo 在 外 推 时 
的 差异 更 小 ， 表 明 加 入 GK 下 的 可 行 性 ， 降 低 了 外 推 的 
不 确定 度 。 


超重 核 “SHN》 的 合成 一 直 是 核 物理 的 前 沿 热点 太 
究 之 一 。 研 究 超重 核 的 途径 主要 是 通过 a 衰变 ， 原 因 在 
T oa 衰变 是 超重 核 的 主要 衰变 模式 。 我 们 基于 潜力 解 解 
集 对 超重 核 CZ > 104) 的 a 衰变 能 进行 了 理论 计算 。 
图 5 为 潜力 解 解 集 对 所 有 超重 核 计算 a 衰变 能 标准 差分 
布 图 。 由 分 布 情况 可 知 此 图 粗略 符合 高 斯 分 布 ， 质 量 最 
优 解 处 于 分 布 中 后 端的 位 置 ， 而 大 部 分 的 潜力 解 集中 在 
质量 最 优 解 的 左 侧 ， 少 部 分 的 潜力 解 在 质量 最 优 解 的 右 
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侧 。 这 表明 在 大 概率 下 港 力 解 所 计算 超重 核 的 a 衰变 能 4 结论 


m 


会 比 质量 最 优 解 更 加 精确 。 另 外 值得 一 提 的 是 ， 原 子 核 


H 


aS 
JU: 


检验 。 


的 壳 演 化 一 直 是 国际 上 的 习 
[ES SEE 
模型 涉及 到 关于 壳 层 结构 的 复杂 价 
(N) =8, 20, 28, 50, 82, N=126, 184 的 壳 层 信息 输入 到 质量 
模型 的 价 核子 计算 之 中 ， 而 实验 中 出 现 的 对 
考虑 到 DZ 模型 的 计算 中 。 本 文 也 对 此 而 
ies rey a a Ia] 


E 点 研究 方向 ， 研 究 表 明和 远离 
B E E 
核子 计算 ， 提 前 将 
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并 对 126 < N < 184 的 a 衰变 能 进行 计算 ， 所 得 结 
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(在 线 彩 图 ) DZ10 模型 以 结合 能 、a 衰变 能 和 GK 关系 


解 。 
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(在 线 彩 图 ) DZ 模型 潜力 解 解 集 计 算 所 有 实验 已 知 的 超 


重 核 CZ > 104) 的 a 衰变 能 标准 差 的 分 布 情况 图 。 黑 
色 箭 头 的 数值 表示 质量 最 优 解 所 计算 的 超重 核 衰 变 能 
标准 差 ， 蓝 色 虚 线 表示 分 布 的 整体 情况 。 
重 核 衰变 能 的 标准 差 ， 纵 轴 为 潜力 解 的 个 数 。 


横 轴 为 所 有 超 


本 文 则 在 通过 添加 严格 的 物理 约束 ， 减 轻 核 质 量 模 
型 的 过 拟 合 程度 ， 从 而 降低 模型 远 距离 外 推 时 巨大 的 不 
确定 度 。 在 多 目标 优化 算法 的 框架 下 , 引入 a 衰变 能 和 局 
域 质量 公式 Garvey-Kelson 关系 , 起 到 大 幅 提 升 oc 衰变 能 
的 精度 ， 增 强 模型 外 推 能 力 的 效果 。 此 外 新 核 素 边界 和 
模型 统计 不 确定 度 通过 多 目标 方法 的 统计 特征 给 出 。 我 
们 从 海量 的 非 劣 解 解 集中 筛选 出 更 具 外 推 能 力 的 解 , ( 称 
为 潜力 解 ) 并 对 潜力 解 做 出 严格 的 统计 分 析 ， 得 出 平衡 
多 个 物理 量 的 全 局 解 比 单一 的 质量 最 优 解 的 预言 能 力 更 
ss 
的 过 拟 合 问题 ， 因 此 提升 了 模型 在 外 推 时 的 可 靠 性 
将 会 对 丰 中 子 核 的 可 靠 预言 提供 一 定理 论 支 持 。 
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Study of multiple physical constraints on the Duflo-Zuker nuclear mass model 


Ye Weihu!, Qian Yibin! 
(1. Department of Applied Physics, Nanjing University of Science and Technology, Nanjing 210094, China) 


Abstract: Nuclear mass plays an important role in both nuclear physics and astrophysics. While the theoretical accu- 
racy of masses has reached a quite astonishing level, the extrapolations among various predictions have been conflicting, 
especially in neutron-rich region. In this paper, we review the main results on the extrapolations of nuclear mass for- 
mulas that we descerned in recent years. By using a rigorous multi-objective optimization on mass models, we take the 
mass differences a decay energy and the Garvey-Kelson relations as multiple physical constraints to reduce the degree 
of overfitting phenomenon, resulting that the predictive power of models was improved in some degree. 

Key words: nuclear mass; a decay energy; multi-objective optimization 
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